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Abstract 

New solutions to the static, spherically symmetric Einstein- Yang-Mills-Higgs equations 
with the Higgs field in the triplet resp. doublet representation are presented. They form 
continuous families parametrized by a = Myj/Mp\ (Mw resp. Mp\ denoting the W-boson 
resp. the Planck mass). The corresponding spacetimes are regular and have spatially com- 
pact sections. A particularly interesting class with the Yang-Mills amplitude being nodeless 
is exhibited and is shown to be linearly stable with respect to spherically symmetric per- 
turbations. For some solutions with nodes of the Yang-Mills amplitude a new stabilization 
phenomenon is found, according to which their unstable modes disappear as a increases 
(for the triplet) or decreases (for the doublet). 

By now a lot is known about classical non-abelian particle-like and black hole solutions in 
Einstein- Yang-Mills (EYM) and in Einstein- Yang-Mills-Higgs (EYMH) theories, we refer to 
for a review and references. The existence of such objects had not been anticipated and their 
study has revealed some important non-perturbative features of the interaction between the 
gravitational and gauge fields. In addition to particle-like solutions where the corresponding 
spacetime is asymptotically flat, a rather different class has been found in the standard SU (2) 
EYMH theory (i.e. without a cosmological constant) with the Higgs field being in the triplet 
representation [2j. They form continuous families parametrized by the gravitational coupling 
strength, a = M\y/Mpi, and the different families are characterized by two integers, the number 
of nodes of the Yang-Mills amplitude, m > 1, resp. of the Higgs field, n > 0. The spacetimes 
corresponding to these solutions have spatially compact sections, i.e. they are topologically 
R x S 3 , representing static non-abelian Einstein-universes. The solutions of Ref. [2j generalize 
those found in EYM theory with a positive cosmological constant A (EYM-A) (Hj. Intuitively 
the self-interaction potential of the Higgs field acts as a "dynamical cosmological constant" 
providing the necessary energy density to support a spatially compact spacetime. The very 
existence of most of these classical solutions is quite surprising and their study is important to 
understand fundamental phenomena in gauge theories coupled to gravitation. 

The important problem of stability of the EYM and EYMH asymptotically flat solutions 
has been studied mostly in linear perturbation theory |U[3]> with the disappointing conclu- 
sion that, with the exception of the gravitating analogues of the 't Hooft-Polyakov magnetic 
monopoles, all of them are unstable. The stability analysis of the compact solutions in the 
EYM-A theory also shows the existence of unstable modes (HI El- Interestingly it has been re- 
cently found that, allowing for a negative cosmological constant, the EYM-A equations admit 
nodeless solutions with asymptotically anti-de Sitter (AdS) spacetimes which have been shown 
to be stable [HI El- 



In this Letter we present globally regular, static and spatially compact solutions in the 
SU(2) EYMH theory with the Higgs field being in the doublet representation, i.e. the non- 
abelian bosonic subsector of the standard model of electroweak interactions. Similarly as for 
the triplet Higgs case [2] , there exist distinct continuous families parametrized by a and indexed 
by the number of nodes of the YM resp. Higgs field, (m,n). We exhibit a new class of compact 
solutions for both the doublet and triplet case where the Yang-Mills amplitude has no nodes. 
These new classes have no analogues in EYM-A theories and their existence is related to an 
interplay of the scalar potential and the gravitational field. They could also be interpreted 
as gravitationally bound (or rather confined) monopole-antimonopole (triplet) or sphaleron- 
antisphaleron (doublet) pairs. 

We have also investigated the stability properties of the compact solutions. The spectrum 
of small fluctuations can be divided with respect to parity into even ("gravitational") and odd 
sectors. A detailed numerical investigation strongly indicates linear stability of these nodeless 
solutions, at least against spherically symmetric perturbations for both the triplet and the 
doublet case in the "gravitational sector" of perturbations which makes their further study 
worthwhile. In the odd parity sector of perturbations an analytic argument shows that for the 
triplet case the number of unstable modes is simply equal to the number of nodes, m, of the 
Yang-Mills field, establishing that the nodeless solutions have no unstable modes at all. In the 
doublet case, however, a suitable adaptation of the result of Ref. leads to the result that 
there is always an unstable mode independently of m. This unstable mode does not couple 
directly to the metric fluctuations and it might just indicate that the static doublet solution 
becomes time-dependent, possibly of oscillatory type since 3-space is compact and there are no 
radiative modes as in an asymptotically flat space-time. 

In addition, we have found yet another interesting phenomenon -"stabilization", which 
occurs in the even parity sector. Our stability analysis has shown that for the triplet case, the 
(m = 1, n = 0) family of compact solutions discovered in [2] has two unstable modes in the 
gravitational sector. Quite surprisingly as a increases the two unstable modes actually merge, 
and for larger values of a no unstable modes are found. This phenomenon can be interpreted so 
that the solutions undergo "gravitational stabilization" since they loose their unstable modes 
for increasing values of the gravitational coupling strength. To the best of our knowledge 
this is the first example that a system gets stabilized when the gravitational coupling gets 
stronger. In fact exploiting the (somewhat unphysical) limit a-*oowe have strong numerical 
evidence that no instabilities occur for arbitrarily strong values of the gravitational coupling. 
Remarkably there is also a family of solutions (m = 2, n = 0) in the doublet case which looses 
its two unstable modes in the gravitational sector of perturbations. This occurs, however, for 
decreasing values of a. The classical solutions we find exist in the regime of strong gravitational 
coupling - at energies close the Planck scale, and therefore they would be mostly relevant for 
semiclassical computations and to gain more understanding of the classical phase space of 
EYMH theories. They seem to be quite distinguished by their stability properties and also by 
the fact that they are present in the standard EYMH theories (i.e. without the introduction of 
new interactions, or negative vacuum energies) making them of quite some interest. 

The most general spherically symmetric line element is 



and the most general spherically symmetric Ansatze for the gauge resp. for the Higgs fields can 
be written as 



A = (A dt + A^dR + cos 9d<p) T 3 + (WT, + KT 2 ) d6 + (WT 2 - KT X ) sin0#, (2a) 





$ = H . 1 + % F T 3 (doublet) , $ = H T 3 (triplet) 



(2b) 
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where T a (a = 1, 2, 3) are the generators of SU(2) and Ai, W, K, H, F are functions of R and t. 
To look for compact solutions we use the "minimal" Ansatz compatible with the field equations, 
corresponding to A = A x = K = F = 0. Following the notations of 0, we define N := r, 
k := N + rz>, where / := df /da := e~ x df /dR, and the static field equations can be written as 

N = -2^-drH 2 , (3a) 

k = [1 - k 2 + 2W 2 - ^(H 2 - a 2 ) 2 - 2C]/r , (3b) 

•• TT W 2 -1 1 dC k — N ■ . , 

W = W — + -— W, 3c 

H 2 oW r 

where [3 = M H /M W , M H being the mass of the Higgs field, and d = 2, C = (W - l) 2 H 2 /2 
resp. d = 1, C = W 2 H 2 for the doublet resp. for the triplet case. It should be noted that 
for the doublet case a is My//y/2Mpi. Additionally there is a constraint equation relating 
{r,N,K,H,W} to the first derivatives of the matter fields. 

Eqs. (j3J) admit a family of local solutions regular at the origin where r = a + 0(a 3 ) and 
N = k = 1 + 0(a 2 ). For the triplet case, the behaviour of the YM resp. Higgs field at the 
origin is W = ±(1 — ba 2 ) + 0(c 4 ) resp. H = aa + 0(cr 3 ) with a, b being free parameters. 
One can impose without restriction W(0) = +1. For the doublet case, however, because of 
the asymmetry of the coupling term, C, between the Higgs and YM fields with respect to 
W — > — W, there are two distinct families of local solutions corresponding to the sign of W(0): 
either W = -(1 - ba 2 ) + 0(a 4 ) and H = aa + 0(a 3 ), or W = +(1 - ba 2 ) + 0(a 4 ) and 
H = a + 0{a 2 ). 

For H (a) = the EYMH eqs. reduce to those of an EYM-A theory with the cosmological 
constant given as A = a 4 /5 2 /4. As found in Ref. [3] for a special set of values of the cosmological 
constant {A = A m , m = 1 . . .} the EYM-A theory admits regular compact solutions. The m = 1 
solution is known analytically jllj : 

r = v / 2sin(cr/\/2) , N = k = W = cos(<r/\/2) , (4) 

with Ai = 3/4. Due to the compactness of the manifold (5 3 in this case) there is an extremum 
of the radial function, r{a e ) = 0, to which we refer as an equator. Following Ref. [2] a natural 
strategy to search for solutions of Eqs. © is to look for solutions bifurcating with those of the 
EYM-A theory. At the bifurcation point the bifurcating solutions are normalizable zero modes 
of the linearized Higgs field Eq. (j2Ji) in the corresponding EYM-A m background solution. In 
the doublet case the equation for the linearized Higgs field amplitude, h, takes the form: 



( W h - l) 2 _Kn\ h _ *b + N h - 

2r^ a 2 J rb 



h=[^r^-^)h-^^h, (5) 



where , iVb • • • etc. stand for the background solution. The bifurcation parameter, a, is to be 
varied to find globally regular solutions of Eq. © . In the analytically known m = 1 background 
(|3J) the bounded zero modes of the Higgs-field, h n , can be explicitly constructed. They are 
characterized by an additional node number, n > 0, and by the values of the bifurcation 
parameter a 2 n = 6/((2n + l)(2n + 5) - 2) , n = 0, 1, . . ., (yielding a 2 = 2, a 2 A = 6/19, . . . ). 
Due to the asymmetric form of C for the doublet, h n is neither symmetric nor antisymmetric 
around the equator, and the solutions of the full system are asymmetric for odd values of the 
Yang-Mills node number. For odd values of m the gauge field amplitude interpolates between 
W = +1 (at the north pole) and W = — 1 (at the south pole) of the compact manifold. 



3 



For even values of to, there are two classes of solutions, both symmetric around the equator 
distinguished by the values of W(0) = ±1 at the origin (referred to as type + resp. — ). We give 
a detailed description of these families in ^2] > here we exhibit a particularly interesting class of 
type W(0) = —1 with m = 2 and n = 0, denoted as (m = 2, n = 0)~ . This solution bifurcates 
with the to = 2 compact solution of the EYM-A theory at a 2fi = 1.494277 with A 2 = 0.364244. 
It exists for a finite range of values of a € [a min , a 2 ,o] with a min « 0.54(9). For decreasing values 
of a, the two zeros of W move towards the equator and they merge there for a* = 0.807469. 
For a < a*, these solutions become nodeless. One can easily show that this cannot happen for 
the (to = 2, n = 0) + solutions. For a — ► a min the spatial sections of the limiting solution will 
no longer be compact (a — * oo). 

For the triplet case, the corresponding analysis has been performed in [2J where various 
one parameter families of solutions of Eqs. I[3jl. indexed by the node numbers to > 1, n > 0, have 
been found. These solutions are symmetric or antisymmetric around the equator depending 
on the parity of to and n. Here we exhibit a new type of family of solutions in the triplet case 
with to = 0, n = 0. This family does not bifurcate with any known EYM-A compact solution 
and it seems to exist for a S [a ,oo) where a 0.87(7). In the strong coupling limit a — > oo, 
Eqs. © reduce to those of an EYMH system with (3 = and with a cosmological constant 
A = 1/4. In this limit the equator, and hence the zero of W, tend to infinity, and the spatial 
sections of the solutions will no longer be compact ( a — > oo and r — * \/2). The limiting metric 
has the asymptotic form ds 2 = a 2 dt 2 — da 2 — 2df2 2 , corresponding to the product of a 2-d 
Minkowski space (in Rindler coordinates) and a 2-sphere. In the limit a — > a the solution 
interpolates between a regular origin and a non-abelian fixed point of Eqs. (j2J) (see Ref. |13j ) 
with non-compact spatial sections. The nodeless solutions for both the triplet and doublet 
cases are displayed on Figs. ^ HI 
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Figure 1: A nodeless triplet solution for a close to a . 

Stability Analysis. We shall restrict our stability analysis to spherically symmetric fluc- 
tuations around the solutions. In many cases instabilities show up already in the spherically 
symmetric sector. The complete stability analysis of the asymptotically AdS EYM-A solutions 
performed in Ref. has demonstrated that solutions stable under spherically symmetric per- 
turbations remain stable under non-spherically symmetric perturbations. The general stability 
analysis is considerably more involved, and it is further complicated for spatially compact solu- 
tions due the existence of an equator. The most general spherically symmetric fluctuations split 
into two sectors according to spatial parity reflections. Making a mode decomposition of the 
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Figure 2: A nodeless doublet solution. 



type f(R,t) = f(R) + Sf(R) e wt , where f(R) denotes a static solution, the linearized EYMH 
Eqs. fall into two sets. In the even parity ("gravitational") sector, the matter perturbations 
(SW, 5H) couple to the gravitational ones {5v, 5 A, 5r). The odd parity ("sphaleronic") sector 
involves only the matter perturbations SAi, 5K and 5F (for the doublet). 

In the even parity sector, there exists a particular gauge - the Schwarzschild gauge - 
defined by 5r = 0, in which the metric perturbations can be eliminated and one obtains a 
closed system of Schrodinger equations for the perturbation amplitudes <I> := (5W , 5H) of the 
form: 



where the 2x2 potential matrix, V, can be expressed explicitly in terms of the functions 
{r, N, e v ,H, W} and their derivatives. For the limiting solutions (a — ► oo) the potential matrix 
is bounded and one can use a generalized version of the Jacobi criterion []3J to detect bound 
states with negative eigenvalues (uj 2 < 0), corresponding to unstable modes. According to this 
criterion, the vanishing of a suitable determinant of the solutions of Eqs. © in (0, +oo) implies 
the existence of at least one bound state with a negative eigenvalue. We have numerically 
integrated © in the background of the limiting (m = 0, n = 0) solution and found that it has 
no unstable mode. 

For solutions with compact spatial sections (i.e. topologically 5 3 ), the potential matrix, 
V, in Eqs. @ is singular due to the presence of an equator, see Refs. El- One can avoid 
such problems by working in a globally regular gauge, e.g. 5X = 0. In this gauge, however, 
one obtains a somewhat unusual eigenvalue problem, which can be nevertheless dealt with 
numerically. Starting at the bifurcation points where the numerical values of the negative 
eigenvalues are known it is possible to determine their dependence on a. We remark that 
the zero mode of the Higgs field at the bifurcation points, becomes an unstable mode for the 
(m = 1, n = 0) triplet solution, comp. Fig. 01 while it goes into the positive spectrum for 
the (m = 2, n = 0)~ doublet solution, comp. Fig. The negative eigenvalues behave quite 
similarly to what has already been observed for certain families in the EYM-A system [7]. It 
is quite remarkable that in the present cases the solutions loose all of their unstable modes. 

For the nodeless triplet solutions, we have undertaken a careful numerical search for 
unstable modes (scanning in the parameter space) and the lack of any evidence for such modes, 
together with the result of the Jacobi criterion in the a — > oo limit strongly suggests that this 
solution has no unstable modes for any a G [«o, oo). 
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Figure 3: Stabilization for the {m = 1, n = 0) triplet solution. The bifurcation point is 
a li0 = 1.732051 

In the odd parity sector, the stability analysis is simplified by the fact that the metric and 
matter perturbations decouple. By a suitable adaptation of the main ideas of Refs. [U Ej we 
have proved the following: a) in the triplet case, solutions with m nodes for the YM amplitude 
have exactly m unstable modes; b)in the doublet case, all solutions have at least one unstable 
mode. 

Therefore we conclude that the nodeless triplet solutions are linearly stable against spher- 
ically symmetric perturbations in both sectors, while the (m = 1, n = 0) triplet solutions sta- 
bilize in the even parity sector but still possess one unstable mode in the odd parity sector. 
The (m = 2, n = 0)~ doublet solutions are stable for a < 1.334 in the even parity sector and 
have at least one unstable mode in the odd parity sector. It is not unlikely that the gravita- 
tional stability persists at the non linear level, and no singularity forms. Then, the instability 
of the static solution indicates merely that it becomes time-dependent (of oscillatory type). 
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Figure 4: Stabilization for the (m = 2, n = 0)~ doublet solution. The bifurcation point is 
a 2i0 = 1.494277 
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